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PREFACE 

The  coordinated  motion  four  bar  linkage  is  a  fairly  common  type  of 
mechanism.   It  has  many  applications.   It  is  the  purpose  of  this  report  to  • 
present  a  method  for  improving  graphical  overlay  designs  by  decreasing  their 
input-output  torque  error.   It  is  also  the  purpose  of  this  report  to  present 
a  method  of  analysis  which  is  more  accurate  and  gives  a  better  indication 
over  the  range  of  operation  of  the  error  between  the  torque  required  to  make 
a  four  bar  mechanism  static,  and  the  torque  produced  by  an.  external  device 
used  to  put  the  mechanism  in  static  equilibrium. 

Indebtedness  is  acknowledged  to  Dr.  John  C.  Lindholm  for  his  suggestion 
for  work  in  this  area.  Also,  I  would  like  to  thank  Dr.  Hugh  Walker  for  his 
advice  on  computer  programming. 
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CHAPTER  I 


INTRODUCTION 


n    Often  it  is  a  necessity  to  design  a  four-bar  linkage  for  which  the  input 
and  output  links  perform  in  some  predetermined  fashion.   The  criteria  for 
excellence  of  performance  is  the  absolute  value  of  the  magnitude  of  error 
with  respect  to  the  desired  input-output  relationship. 

The  problem  considered  is  one  for  which  the  desired  relationship,  which 
may  be  determined  from  external  torque  requirements,  involves  producing 
coordinated  input-output  crank  motions. 

A  graphical  method  for  design  of  a  four-bar  mechanism  producing  angular 
displacements  according  to  an  input-output  schedule  is  available  and  quite 
useful.   Seeking  to  reproduce  the  graphical  method  on  the  computer  led  to 
a  variation  of  parameter  technique  which  aids  the  linkage  designer  multifold. 

The  Graphical  Overlay  Method 
In  order  to  design  a  four-bar  mechanism  which  will  produce  a  specified 
output  angle  for  a  specified  input  angle  a  simple  procedure  must  be  followed i 

1.  Develop  a  schedule  of  positions  for  the  specified  input  angles  and 
the  related  output  angles . 

2.  On  tracing  paper  make  a  layout  (Figure  1)  showing  the  successive 
positions  of  the  first  crank.   Choose  any  desirable  length  for  this 
crank.   Assume  a  length  for  the  connecting  rod,  and  draw  the  family 
of  circular  arcs  of  this  radius  with  their  centers  at  the  successive 
crankpin  positions  (Figure  2). 
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Figure  1,  Layout  of  Successive  Positions 
for  the  First  Crank 
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Figure  2,  Successive  Positions  for  Crank  One 
and  Loci  of  Positions  for  Extreme  End  of 
Connecting  Rod 
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3.  Make  a  second  layout  (Figure  3)  showing  the  successive  positions  of 
the  second  crank  and  a  series  of  possible  lengths  for  this  crank. 

4.  Fit  the  first  layout  over  the  second,  as  shown  in  Figure  4,  trying 
to  make  the  circular  arcs  of  the  first  pass,  in  proper  order, 
through  one  of  the  series  of  constant  radius  crankpin  positions  of 
the  second.   It  may  not  be  possible  to  obtain  a  completely  satis- 
factory fit  with  the  first  choice  for  a  connecting  rod.   If  this  is 
the  case,  the  first  layout  must  be  redrawn  with  a  different  choice  of 
length  for  the  connecting  rod,  [l]*  pp.  45-46. 

It  is  generally  accepted  that  the  accuracy  of  such  a  method  is  in  the 
order  of  +  .5  degrees.  The  only  method  for  improving  accuracy  is  to  construct 
larger  overlays  containing  more  intermediate  positions.  Calculation  of  the 
error  with  respect  to  the  desired  relationship  is  extremely  tedious  and  errors 
in  calculation  become  cumulative.  The  end  result  is  a  calculated  value  of 
error  which  may  have  very  little  significance.  Since  nothing  is  known  about 
the  error  between  the  positions  chosen  for  the  overlays,  these  calculations 
may  not  even  indicate  the  true  value  of  maximum  error. 

The  graphical  method  offers  no  control  over  link  lengths  and  it  often 
becomes  awkward  trying  to  satisfy  as  many  intermediate  positions  as  possible 
over  a  specified  range. 

Two  good  features  of  the  graphical  method  are  that  it  contains  no  com- 
plicated theory  and  it  is  easy  to  construct  the  overlays  which  will  produce 
fairly  good  results.  Essentially,  the  graphical  method  varies  three  parameters 
simultaneously.  They  are:  link  one,  link  four,  and  the  initial  position. 


*See  Selected  Bibliography 
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Figure  3,  Second  Layout  with  Positions  of 
Crank  Two  Corresponding  to  Required 
Positions 
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Figure  4,  Overlay  of  Layouts  to-  Determine 
Length  of  Crank  Two  and  Link  One 

A.  Crank  One  (Link  Two) 

B.  Crank  Two  (Link  Pour) 

C.  Connecting  Rod  (Link  Three) 

D.  Distance  Between  Centers  (Link  One) 


.. 


The  designer' 8  eye  is  used  as  the  feedback  mechanism  to  produce  least  error. 
Should  the  design  appear  to  be  unsatisfactory  the  designer  changes  the  value 
of  link  three  and  repeats  the  graphical  overlay  procedures. 


The  Variation  of  Parameter  Technique 
There  are  five  important  parameters  assaciated  with  working  a  coordinated 

motion  four  bar  linkage  problem.  They  are  the  lengths  of  the  links  one 

■ 

through  four  and  the  initial  position  of  the  mechanism  (see  Figure  B-l  in 
Appendix  B) . 

After  having  drawn  several  graphical  solutions,  it  was  observed  that  the 
chances  for  obtaining  a  completely  satisfactory  mechanism  on  the  first  trial 
were  roughly  50  percent.  No  positive  evidence  was  given  by  the  graphical 
process  to  support  a  claim  that  a  better  solution  could  be  obtained  only  by 
increasing  or  decreasing  the  value  of  some  parameter. 

It  seems  heuristically  evident  that  if  a  solution  which  is  close  to 
having  minimum  error  is  obtained,  a  parameter  at  a  time  may  be  varied  over 
a  range  containing  the  base  solution;  and  if  the  parameter  variation  which 
produces  the  least  absolute  value  of  maximum  error  is  retained  as  a  new  base,  a 
minimum  error  will  be  approached.-  While  it  is  true  that  the  best  path  of  var- 
iation may  not  have  been  chosen  (if  the  error  is  dependent  on  path)  a  better 
than  initial  mechanism  is  always  produced. 

If  the  absolute  value  of  error  versus  a  particular  parameter  (link  length 
or  initial  position)  is  graphed  continuously  for  a  range  of  values  containing 
the  base  design,  the  validity  of  the  argument  can  be  easily  seen.  Figure  5 
shows  a  hypothetical  case  for  which  the  absolute  value  of  error  is  plotted  against 
a  range  of  a  particular  link. 

In  this  example  the  error  is  calculated  for  nine  link  lengths.  One 
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Figure  5 ,  Error  for  a  Particular  Link  With 
All  Other  Parameters  Held  Constant 

A.  Value  of  Link  Chosen  as  a  Minimum 

B.  Base  Design  Link  Length 
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of  these  calculations  must  be  the  base  design.  The  minimum  calculated  value 

of  maximum  absolute  error  is  for  link  value  A.   It  should  be  noted  that  the 

i 

error  at  length  A  may  not  be  the  absolute  minimum  in  the  range. 

As  the  method  directs,  some  other  parameter  is  optimized  in  the  same  fash- 
ion. If  the  maximum  absolute  value  of  error  versus  the  same  link  for  which  the 
error  was  first  minimized  is  graphed,  it  is  noted  that  the  new  absolute  value 
of  error  for  link  length  A  is  either  as  small  of  smaller  than  the  minimum  error 
due  to  the  first  parameter  optimization.  ... 

The  steps  for  solving  the  problem  after  the  desired  relationship  between 
input  and  output  crank  is  established  are  as  follows:     »'..'".. 

1.  Develop  a  graphical  solution. 

2.  Choose  a  small  range  of  values,  containing  the  base  design  values, 

in  which  some  or  all  of  the  parameters  are  to  be  varied.  The  design- 
er may  choose  not  to  vary  some  parameter. 

t 

3.  Vary  a  single  parameter  at  a  time  over  its  entire  range.  The  other 
parameters  are  constants  during  this  operation. 

4.  Retain  the  parameter  value  which  renders  the  least  absolute  value 
of  maximum  error  as  a  new  base. 

5.  Vary  any  other  parameter  as  prescribed  in  3  and  4.* 

6.  Repeat  the  process  until  the  error  can  no  longer  be  reduced  or  a 
satisfactory  result  is  obtained. 

The  most  evident  advantage  of  the  method  is  the  designer  has  full 
control  over  as  many  parameters  as  he  wishes  to  control.  If  done  on  the 
computer,  the  process  can  be  made  to  compute  the  error  at  many  intermediate 
positions,  with  tremendous  speed,  and  great  accuracy.  Thus,  more  can  be 
learned  about  the  original  graphical  design  and  its  subsequent  improvements. 
The  program  for  doing  this  is  discussed  in  Chapter  IV. 
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CHAPTER  II 

DEFINITION  OF  THE  PROBLEM 

To  further  define  the  type  of  problem  considered,  It  would  be  wise  at 
tt.is  time  to  state  some  assumptions.  They  are: 

1.  All  links  are  rigid. 

2.  All  links  are  weightless. 

3.  The  linkage  is  static. 

4.  All  solutions  are  limited  (due  to  complexity  of  geometry)  to  the 
region  above  the  infinite  line  through  link  one's  center line. 

5.  Link  two  plus  link  three  plus  link  four  must  be  greater  than  link  one. 

6.  Only  pure  torques  are  applied  to  the  mechanism.  They  are  on  link 
two  and  link  four. 

7.  The  pins  are  frictionless  in  the  analysis. 

8.  Neither  of  the  cranks  can  reverse  direction  over  the  prescribed  range. 

9.  Both  cranks  travel  either  counterclockwise  or  clockwise. 

With  these  assumptions  it  is  possible  to  analyze  the" four  bar  extensively. 
An  analysis  is  presented  in  Appendix  A. 

If  there  is  an  input-output  torque  relationship  which  must  be  (maintained, 
one  torque  function  may  be  applied  to  either  the  input  or  output  crank  and  the 
parameters  may  be  adjusted  so  that  the  balancing  torque  on  the  other  crank 
approaches  the  second  torque  function  over  a  specified  range. 


: 
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CHAPTER  III 

A  SPECIFIC  PROBLEM 

To  further  illustrate  the  variation  of  parameter  process,  consider  the 
fo  lowing  problem,  which  is  problem  E3.5  in  Hall's  text  [l]. 

Design  a  spring-linkage  counterbalancing  system  of  the  type  shown  in 
Figure  6.  The  weight,  W,  is  50  lb  and  the  pulley  radius  is  12  in."  It  is 
required  that  the  weight  be  counterbalanced  over  a  90-deg.  range  of  pulley 
rotation,  e.g.,  within  this  range  the  weight  should  remain  in  any  position  in 
which  it  is  placed.  Neglecting  friction,  this  would  mean  that  the  force 
exerted  by  link  3  on  link  4  should  have  a  moment  about  the  pivot  just  equal 
to  the  moment  WR.   However,  it  is  anticipated  that  there  will  be  enough  fric- 
tion in  the  joints  so  that  it  will  be  necessary  to  balance  the  moment  WR  only 
within  +  3%. 

The  spring  is  assumed  to  be  linear  and  is  to  be  assembled  so  that  the 
spring  torque  (exerted  on  link  2)  will  range  from  T.    to  2T.  over  the  counter- 
balancing range.  Design  for  a  90-deg.  range  of  motion  for  link  2. 

Solving  for  the  spring  constant  (k) : 

k  -  (2Tt  -  T^Q  -  T^AO  -  ^/OT/2)  -  2Tt/lT. 

The  torque  in  the  spring  is : 

T  -  T  +  Ok. 

The  potential  energy  lost  due  to  a  downward  displacement  of  the  weight  is 

stored  in  the  helical  spring.  By  using  an  energy  balance,  the  initial  torque, 

T^,  and  the  functional  relationship  between  input  and  output  displacement 

angles  is  established. 

The  energy  balance  is: 


KT.    +  0k)d9  -   )WRd0. 

Solving  for  T,: 

T    -  100/3  foot  pounds. 
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Figure  '6,  General  Picture  of  the 
Example  Problem   [l] 
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Solving  the  energy  equation  for  some  intermediate  position  6  and  0  : 


! 


°»  ,« 


■  • 

! 


n 


(T^  6k)d0  «=  JWRd0; 


b 

TiGn  +  kO^/2  -  WR0n. 
From  this  equation  a  table  for  input  and  output  angular  displacements  was 
developed,  Table  1.  A  graphical  solution  was  obtained.   It  is  displayed  as 

Figure  7.   Note  that  step  one  has  just  been  completed,  [pf; 

,  ■ 

The  length  of  link  four  and  the  initial  position  were  the  parameters 

chosen  to  vary.   The  values  and  ranges  for  these  parameters  were: 

j  ■. 
Link  four  -  2.06  +  .06  units; 

output  angle  =  .6455  +  .008  radians. 

This  completed  step  two. 

■ 

Starting  with  the  smallest  value,  link  four  was  positively  incremented 
until  step  four  was  complete.   The  output  angle  was  incremented  similarly. 
For  this  problem  one  cycle  of  these  two  parameter  variations  was  sufficient   • , I 
to  illustrate  the  process. 


The  program  for  doing  this  is  in  Appendix  C.   An  increment  size  of  .020 
units  was  chosen  for  link  four.   An  increment  size  of  .002  radians  was  chosen 
for  the  initial  output  angle.   The  computer  printout  for  the  problem  is  shown 
in  Table  2. 

As  the  computer  incremented  the  value  of  link  four  through  its  range,  it 
n-ut  only  calculated  the  error  for  the  graphical  design,  but  also  computed  the 
error  for  a  link  four  length  that  had  less  error  than  the  graphical  design. 

The  best  linkage  in  the  range  was  when  link  four  equaled  2.02  units. 
The  maximum  error  was  computed  to  be  7.857  percent  for  this  link  length.   The 
graphical  design  had  a  maximum  error  of  10.99  percent.  This  was  an  improve- 
ment of  3.033  percent  due  to  only  one  parameter  variation  process.   Since  the 


TABLE  1 
SCHEDULE  OF  POSITIONS  FOR  THE  EXAMPLE  PROBLEM 
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Position 

©   (degrees) 

o   (radians) 

jJ     (degrees) 

4>    (radians) 

0 

0.00 

0.0 

0.0 

0. 

',     1 

11.25 

tr/16 

6.99 

0.1219 

2 

22.50 

tr/8 

16.90 

0.2945 

3 

33.75 

3*716 

26.65 

0.4660 

.     4 

45.00 

TT/4 

37.45 

0.6540 

5 

56.25 

5*716 

.     49.20 

0.8590 

6 

67.50 

3*78 

61.90 

1.0800 

7 

78.75 

7  IT/16 

75.50 

1.3180 

8 

90.00 

tr/2 

90.00 

1.5710 

■  ■ 
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Figure  7,  A  Graphical  Solution  for 
the  Example  Problem 


A. 
B. 
C. 
D. 

E. 
F. 


Link  One  = 
Link  Two  ** 
Link  Three 
Link  Four  « 


3.19 
2.50 
=  A. 085 
■  2.06 


Initial  Angle  (A  +  B)  - 
Initial  Output  Angle  -  , 


2.9693  radian 
6455  radian 


\ 
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COMPUTER  PRINTOUT  FOR  THE  EXAMPLE  PROBLEM 


LINK  1=  O.J189Q990fZ  6l     LU<K     .>  =   0.25000000E  01 

LINK  3=  O.4O8S00OUF  01    LINK  4  =   0.20000000E  01 

INITIAL  INPUT  ANGLE=   0  •  29  1  30  350E  01 

INITIAL  OUTPUT  ANGLE=   0.64550000E  00 

MAXIMUM  ERRGR=   0.97238420E-0 1 


LINK  1=  G.31899990E  01     LINK  2=   0.2SOOO00OE  01 

LINK  3=  0.40850000E  01    LINK  4  =   0.20199990E  01 

INITIAL  INPUT  ANGLE=   0.292754S0E  01 

INITIAL  OUTPUT  ANGLE=   0.645SOOOOE  00 

MAXIMUM  ERRGR=   0. 78S7S780E-0 1 

NOT  LEAST  ERROR  WHEN  LINK  4  =   0.20399990E  01 

LINK  1=  0.31899990E  01    LINK  2=   O.2SOOO0OOE  01 

LINK  3=  G.40850000E  01    LINK  4  =   0.203999<J0E  01 

MAXIMUM  ERHOR=   0. 05260800E-0 1 


NOT  LEAST  ERROR  WHEN  LINK  A=   0.20599990E  01 
LINK  1=   0.31899990E  01    LINK  2=   C25000G00E  0 V 
LINK  3=   O.A0850000E  01    LINK  4  =   0.20599990E  01 
I.      MAXIMUM  ERRORS   0.10992690E  00 


NOT  LEAST  ERROR  WHEN  LINK  4=   0 . 20799990E . 0 1 
LINK  1=   0.31899990E  01    LINK  2=   O.2bOO0O00E  01 
LINK  3=   0.A0850000E  01    LINK  4=   0.20799990E  01 
MAXIMUM  ERROH=   0.13346610E  00 


NOT  LEAST  CRRQH  WHFlN  LINK  4=   C.20999990E  01 
LINK  1=   0.31899990E  01    LINK  2=   0.25000000E  01 
LINK  3=   0.40850000E  01    LINK  4=   0.20999990E  01 
MAXIMUM  ERROR=   0.1SS91940E  00 


NOT  LEAST  ERROR  WHEN  LINK  4=   0.21199990E  01 
LINK  1=   0.31899990E  01    LINK  2a   0.25000000E 
LINK  3a   0.40850000E  01    LINK  4a  .  0.2U99990E 
MAXIMUM  ERRCRa   0.17733370E  00 


'■  .• 


•  •  ■'■ 
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LEAST  ERROR  DUE  TO  LINK  4 

LINK  1=   0.31899990E  01    LINK  2-        0.25000000E  0 

LINK  3=   0.40850000E  01     LINK  4  =   0.20199990E  0 

INITIAL  INPUT  ANGLE=   0.2927b4'jOI:  01 

INITIAL  OUTPUT  ANGLE=   0 .  6  V550000E  00 

MAXIMUM  ERRCR=   0  •  7tl57&780E-0  I 


LINK  I  =  ,   0.318999O0E-  01     LINK  2  =   0.2'jOOOOOOE  0 

LINK  3=   0.40850000E  01    LINK  4  =   0.20199990E  0 

INITIAL  INPUT  ANGLE=   0.29394910E  01 

INITIAL  OUTPUT' ANGLE=   0.63749990E  00 

MAXI^Uf  ERRCR=   0.75980120E-01  .  , 


■ 


LINK  1=  0.31899990E  01     LINK  2=   0.2'iOOOOOOE 

LINK  3=  0.40850000E  01    LINK  4=   0.20199990E 

INITIAL  INPUT  ANGLE=   0.29364900E  01 

INITIAL  OUTPUT  ANGLE=   0.63949990E  00 

MAXIMUM  ERROR=   0. 73909870E-0 1 


01 
01 


NOT  LEAST  ERROR  WHEN  OUTPUT  ANGLE: 
LINK  1=   0.31899990E  01     LINK  2= 
LINK  3=   0.408S0000E  01    LINK  4= 
MAXIMUM  ERR0R=   0. 75470860E-0 1 


0.64149990E  00 
0.25000COOE .01 
0.20199990E  01* 


NOT  LEAST  FRROR  WHEN  OUTPUT  ANGLE=   0.64349O90E  00 


LINK  1=  0.31899990E  01  LINK  2= 
LINK  3=  0.40850000E  01  LINK  4= 
MAXIMUM  ERRCP=   0 • 7 7025890E-0 1 


0..2SOOO0OOE  01 
0.20199Q90E  01 


1 


NOT  LEAST  ERROR  WHEN  OUTPUT  ANGLE=  0.64549<J90E  00 

LINK  1=   0.310999O0E  01    LINK  2=  0.25000000E  01 

LINK  3=   0.40850000E  01    LINK  4=  0.Z0199990E  01 
MAXIMUM  ERROR=   0  .76575780E-0 1 


NOT  LEAST  ERROR  WHEN  OUTPUT  ANGLE' 
LINK  1=   0.31899990E  01     LINK  2= 
LINK  3=   0.40850000E  01    LINK  4= 
MAXIMUM  ERRCR=   0  •  80 1 1 8230E-0 1 


0.64749990E  00 
0.25000000E  01 
0.20199990E  01 


NOT  LEAST  ERROR  WHEN  OUTPUT  ANGLE*   0.64949990E  00 

LINK  1=   0.318999O0E  01    LINK  2s  0.25000000E  01 

LINK  3»   0.408S0000E  01    LINK  4*  0.20199990E  01 
MAXIMUM  ERRGRs   0  .0  1  69 7400E-O I 
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NOT  LEAST  EHROR  WHEN  OUTPUT 
LINK  1=   0.318999008  01    LI 
LINK  3=   0.A0850000E  01    LI 
MAXIMUM  ERPGW=   0.633I7330E- 


: 


ANGLE*   0.6bl49990E  .00- 

NK  2  =  0.25000000E  01 

NK  "4  =  0.20199Q90E  01 
01 


■":  • 

- 


NOT  LEAST  EWrtOR  WHEN  OUTPUT 
LINK  1=   0.31899990E- 01    LI 
LINK  3=   O.408S00O0E  01    LI 
MAXIMUM  ERROR=   0.84927490E- 


ANGLE=   0.65349990E  00 

NK  2s  0.25OQO0O0F  01 

NK  4  =  0.20199990E  01 
01 


LEAST  ERROR  DUE  TO  OUTPUT  AN 
LINK  1=  0.31899990E  01  LI 
LINK  3=  0.408SOOOOE  01  LI 
INITIAL  INPUT  ANGLE=  0.2936 
INITIAL  OUTPUT  ANGLE=  0.639 
MAXIMUM  ERR0R=   0.73909870E- 


CLE 

NK  2=   O.2bO0000OE  01 

NK  4=   0.20199990E  01 

4900E  01 

49990E  00  .  « 

01 


<'..  ' 


. 

■   ■ 
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maximum  error  was  less  when  link  four  equaled  2.02,  this  value  was  retained 
as  the  new  length  for  link  four. 

Next,  the  output  angle  was  varied  as  prescribed.   The  least  maximum  error 
for  this  variation  process  was  found  to  be  7.391  percent  when  the  initial 
output  angle  equaled  .6395  radians.  This  represents  an  overall  reduction 
of  the  maximum  error  of  3.599  percent  and  ".466  percent  reduction  of  the  pre- 
vious least  maximum  error. 

Further  similar  parameter  variations  will  decrease  the  error  even  more. 
The  variation  of  these  two  parameters  was  executed  in  48.78  seconds  on  the  IBM 
360/50  computer. 

This  program  produced  16  parameter  combinations.  For  each  parameter 
combination  the  error  was  computed  at  every  two  degrees  over  a  90  degree 

y 

* 

displacement  of  the  output  angle. 

■ 
.... 


1 


/■  •    ■ 

CHAPTER  IV 

THE  PROGRAM  AND  ITS  USE 


It  is  assumed  that  the  designer  has  at  least  a  cursory  knowledge  of  Fortran 
programming  techniques.  This  particular  program  was  compiled  and  executed  on 
the  IBM  360/50  computer.  The  IBM  360/50  accepts  cards  punched  on, an  IBM  29 
Card  Punch. 

The  program  is  essentially  divided  into  four  sections.  They  are: 

■ 

1.  a  data  section, 

2.  a  variation  of  link  four  and  a  variation  of  initial  output  angle 
section, 

3.  an  input  torque  subroutine  section,  and 

-1 
A.   two  sin  (argument)  functions.'  .-..:■. 

The  data  section  is  a  group  of  cards  immediately  following  the  format 

cards.  The  symbols  and  their  interpretations  are  as  follows: 

1.  XLINC  is  the  value  of  the  increment  size  for  the  link  length.  It 
may  be  any  value  such  that  when  it  is  multiplied  by  some  integer 
and  is  added  to  CXL4,  the  base  design  value  for  the  linkage  is 
obtained. 

2.  AINC  is  the  value  of  the  increment  size  for  the  initial  output  angl 
It  is  measured  in  radians.   It  may  have  any  value  such  that  when 

it  is  multiplied  by  some  integer  and  is  added  to  AANG,  BANG  is 
obtained. 

3.  AANG  is  a  value  of  the  initial  output  angle  which  is  smaller  than 
the  graphical  design  value.  It  is  measured  in  radians.  It  must 
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have  a  value  such  that  an  integer  times  AINC  plus  AANG  equals  BANG. 

4.  BANG  is  the  graphical  design  value  for  the  initial  output  angle. 
It  is  measured  in  radians. 

5.  XL1  is  link  one. 

6.  XL2  is  link  two. 

7.  XL3  is  link  three. 

■     ■  i-  naj; 

8.  CXL4  is  a  constant  value  for   link  four.     It  is  smaller  than  the 

;  i,  < '    ■  iion    ■ 

graphical  design  value  for  link  four  and  has  a  value  .such  that 

when  it  is  added  to  the  product  of  an  integer  and  XLINC,  the 

graphical  value  of  link  four  is  obtained. 

■  "'Uiah 

9.  SK  is  the  spring  constant.  Omit  this  card  if  there  is  no  spring 

on  the  input  crank. 
10.  TI  is  the  initial  torque  on  the  spring.  Omit  this  card  if  there 
is  no  initial  torque. 

■  1  * 

The  variation  of  link  four  and  initial  output  angle  section  follows  the 
data  section.  The  output  angle  is  incremented  in  2T7/180  radian  steps  over 
afT/2  radian  displacement.  This  may  easily  be  changed. 

In  order  to  change  the  increment  size  of  the  output  angle  simply  change 
the  card  in  this  section  incrementing  the  output  angle  from  XOUTA  =  XOUTA  :  T 
+  2.*3. 141592/180.  to  XOUTA  -  XOUTA  +  desired  increment  size  in  radians.  Th« 
output  angle  displacement  divided  by  the  desired  increment  size  must  be  an 
integer  (N) .  The  DO  loop  index  must  be  changed  to  I  -  1,  N  +  1. 

The  torque  the  designer  is  trying  to  simulate  is  the  torque  developed  in 
a  spiral  spring  (TSP).  The  designer  may  wish  to  use  another  device  on  the  in- 
put crank.  This  will  necessitate  the  replacement  of  the  TSP  card  by  an  ap- 
propriate function.  Also,  the  error  (A5)  card  must  be  repunched  to  conform 

■ 
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C 


;■   •  '  ■l).-'n  '(-!  . 


:„,,„ 

to  the  new  expression  for  error.  {»|1(; 

If  it  is  desired  to  vary  more  parameters,  the  inserted  cards  follow  a, 
pattern  similar  to  the  pattern  established  by  the  first  two  parameter  varia- 
tion processes.  •   *^ 

The  torque  subroutine  calculates  an  input  torque  for  a  constant  output  ;u_ 
torque.   The  output  torque  may  be  easily  changed  from  a  constant  torque  to  a 
torque  with  a  value  dependent  on  the  output  angle  by  replacing  the  card  (  I|Vrt 

■ 

marked  TOUT  =50.  to  TOUT  -  f(XOUTA). 

Since  there  is  a  continuous  need  to  solve  for  angles  in  radians,  funcr_... 
tions  ANGLE  and  GANGL  were  developed,  ANGLE  may  be  used  for  angles  from 
0  to  90  degrees.   GANGL  is  used  for  angles  from  90  to  180  degrees.   To  use    .  . 
either  of  these  functions  simply  submit  the  sine  of  an  angle  to  the  appropri- 
ate function  and  the  value  of  the  angle  will  be  calculated  and  returned  to 
the  main  program. 


'  ".\ 


,  \  . 
•  ■ 


■     ', 


CHAPTER  V 
SUMMARY  AND  CONCLUSIONS 

The  steps  involved  in  using  the  variation  of  parameter  technique  are 
outlined  clearly  in  Chapter  I.  The  steps  are  simple  to  follow  and  are  aided 
greatly  by  the  computer  program  listed  in  Appendix  C.  The  use  of  .this  program 
is  described  in  Chapter  IV. 

It  was  demonstrated  that  the  error  can  be  significantly  decreased  by 
applying  the  variation  of  parameter  technique.  In  the  example  problem  the 
graphical  error  was  10.99  percent.  This  error  was  reduced  to  7.391  percent 
after  applying  the  variation  of  parameter  technique  first,  to  link  four,  and 
second,  to  the  initial  output  angle.  Further  applications  of  the  technique 
will  reduce  the  error  even  more. 

These  two  applications  were  executed  in  48.78  seconds  on  the  IBM  360/50  . 
computer.  This  represents  46  calculations  of  the  error  for  each  of  16  para- 
meter combinations.  Because  the  error  was  calculated  in  many  intermediate 
positions,  the  maximum  error  the  linkage  produced  was  accurately  determined. 

The  number  of  intermediate  positions  can  be  increased  and  the  size  of 
parameter  incrementation  can  be  decreased  to  produce  an  even  better  analysis. 

From  the  results  of  the  example  problem  it  is  concluded  that  the  varia- 
tion  of  parameter  technique  can  be  employed  t6  reduce  error  and  yield  a  design 
which  is  more  satisfactory  than  the  design  produced  by  graphical  method. 
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Appendix  A 
Equilibrium  Equations 
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Figure  A-.l,  Torque  and  Force  Components 
on  a  Separated  Four  Bar  Mechanism 

Separating  the  mechanism  at  its  pins  as  in  Figure  A-l,  the  force 


*  •  *    . 

..  ■:• 


balance  for  each  link  is: 

0 


F   e-.P   •  0; 
14x   34x 


F7/  -  Fi/  "  0; 
34y   14y 

F43x-P23x"0! 

'«Jy-  P23y"  0i 

F32»"  f12x"  °- 


., 

• 

■ 

'" 

• 

: 

'        ■ 
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•• 


■ 


The  force  balance  for  each  pin  isj 
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•0 


F34x"  ?43xJ 
F34jTF43y; 
F23y"  F32y; 
F23x"F32x* 


■  ■  ■ 

2ft  \ 
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Pigure  A-2,  Forces  and  Torque  on 
Links  Three  and  Four 

Performing  a  balance  of  torques  about  point  "A"  for  Figure  A-2  obtains: 


Tout  ut"  F23x^Link  A><8tn(0utPut  AnSle)  -  (Link  3)(aln(K)Jr 
it      /jLink  4)(co»  Output  Angle)  +  (Link,3)(coa(lQ7  -  0. 
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Figure  A-3,  Forces  and  Torque  on  Link  4 
Performing  a  balance  of  torques  about  point  "A"  for  Figure  A-3,  and 


realizing  F3^x  -  F23x  and  F34y  -  F23y  obtains: 


c 


• 


' 


Toutput  -  p23x(Llnk  4)(sin(0utput  Angle) 


•F23y(Link  4) (cos (Output  Angle)  ♦  0. 


Let: 


■ 


BJ 


• 


AA  -  (Link  4)(sin(Output  Angle)), 
.  BB  -  (Link  3)(sin(K)), 

CC  -  (Link  4) (cos (Output  Angle)); 

DD  -  (Link  3)(cos(K)). 
Equation  l_j^J   becomes: 

F23x(AA  -  BB)  +  F23y(CC  +  DD)  -  Toutput, 
Equation  [_2-J  becomes : 

F23*^>  +F23y(CG)  -Toutput. 
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Solving  for  F23x  and  ?23y: 


'F 


23x 


Loutput 


Loutput 


CC  +  DD 
CC 


m     output 


(DD) 


F 


23y 


AA  -  BB 

CC  +  DD 

(BB)(CC)  +  (AA)(DD) 

AA 

CC 

• 

. 

AA  -  BB 

T 
.   output 

• 

■ 

AA 

•  ^output 

T   ■   „  (BB) 
■  output  x   '         , 

■ 

o& 


(.. . 


\J 


AA  -  BB 

AA 


CC  +  DD 
CC 


(BB)(CC)  +  (AA)(DD) 


OJ. 


•  ' 


■  ■ 


• 


/ 


■ 


\    F32y 


F32x 


4 

,  -  .- 

'  •  ■  ■  ■  . 

■  ■ ,  •• 

'  '  '-I 


j 


input 


Figure  A-4,  Forces  and  Torque  on  Link  2 
Performing  a  balance  of  torques  about  point  "B"  for  Figure  A-4,  and 
realising  F32x  -  F23x  *nd  F32y  "  F23y  yields: 

Tinput  -  F23jc(Link  2)(sin(A  +  B))  -  F23y  (cos(A  +  B)). 
It  is  possible  at  this  point  to  solve  for  Tinput  for  any  Tout?ut. 
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Appendix  B 

Geometry 
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Link  3 
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Figure  B-l,  Geometry  of  a  Four  Bar  Linkage 
If  all  the  links  and  G  +  H  are  known: 

Link  5  -  sqrt/(Link  *)2+  (Link  l)2 -2 (Link  l)(Link  4)cos(G  +  1$ 
A  -  s in"1/!  in (G  +  H)(Link  4) /Link  5/ 

C-H)  -  cos'^^Link  5)2-(Link  2)2-(Link  3)2_7/  /,-2(Link  3)(Link 
B  -  8in*1/7sin(C-«))  (Link  3)/Link  5 J 

■ 

E  -  8in-l/.~sin  (C-H))  (Link  2) /Link  5_7 

F  -  8  in"  ^"s  in  (G-fH)(Link  1)/Link  5_7 

Link  6  -  sqrt  /""(Link  l)2+(Link  2)2-2(Link  l)(Link  2)cos(A4B)_7 

G  -  ain-^'sin   (E-fF)    (Link  3)/Link  6J 

H  -  s  in"  1Lb in  (A4B)  (Link  2) /Link  6 J'" 

C  -  ain-l/.-Bin  (H)  (Link  1)/Link  2 J 
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D  -  sln-l/f8111  (E-ff)    (Ltnk  4)/Llnk  hj 
K  -  180°   -   (G-W)   -   (E4F) 
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Appendix  C 


Computer  Program 
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9000  F0WHAT(BH  LINK  1  =  ,  ET  1  to  .  fl  .  2X  ,  8H  LINK  2  =  ,  E  I  6  .  Q..'./  ,  «H  LINK  J=.E16. 
12X,bH  LINK  4=,Elto.8,/,21H  INITIAL  IM'UT  ANGLE=  ,  E  1  6  •  0  ,  /  . 
222H  INITIAL  OUTPUT  ANGLE=  ,F.  I  6.  8  ,  /  ,  I  5H  MAXIMUM  ERROR?  ,  E  16  •  8  ,// 

9002  FQRMATIbH  LINK  1  =  ,  E  1  b  .  8  ,  2  X  ,  8H  LINK  2=  .  E  I  0  .  8  .  /  ,  8H  LINK  3*.Elb. 
!2X,fiH  LINK  4=,Elb.8,/«  1  5H  MAXIMUM  LKWOR  =  ,  E  1  6  .  6  ,  //  ) 

9003  FORMAT (26H  LEAST  ERROR  DUE  TO  LINK  A) 

9004  FURMATU7H  VO 1 0  SET  LINK  4=,ri6.8.//)  •  '    c 

9005  FORMAT (29H  NOT  LEAST  ERROR  WHEN  LINK  4=,Elb.8) 
.9007  F0RMATI35H  NOT  LEAST  ERROR  WHEN  OUTPUT  ANGLE* . E 1 6 .6 ) 

9008  F0RMATI23H  VOID  SET  OUTPUT  ANGLE= ,E 16. 8t // ) 

9009  FORMAT (32H  LEAST-  ERROR  DUE  TO  OUTPUT  ANGLE) 


XLINC=.020 
AINC=.002 


.   .  .   .  . 

.     ..."    « 


■ 
•   : 


AHA=0. 

4ANG=.G375 

FIANG=.6A55 


ii  ' 

XL1=3.19 


XL2=2.50 


XL3=4.0U5  •        ''■,  :■;•',? 

CXLA=2.00- 


XL4=CXL4 


1  SK=200./(3.*3. 141592)  ' 


' 


• 


- 


2  TI=100./3. 
RN0=2000o 
UNO2=20O0.  ..     \  .;  .  .    • 

100  XOUTA=DANG 

101  XFRR=0. 
COU=COU+l. 

CALL  TORO  <XL1.XL2.XL3,XL4,X0UTA,TIN.CSB,CSA,SNAB,CSAB.E,F. 
1SNA,SNH.XKR,F23X,F23Y,A.U)  . 

17  ZERO=A*B  . 

ZOUTA=XOUTA 

00  300  1=1, A6 

CALL  TORQ  ( XL  I . XL2 . XL3, XL4 . XOUTA , Tl N ,CSB , CSA . SNAO  ,CSAB ,E ,F , 
y>        1SNA,SNH,XKR,F23X,F23Y,A,B) 

TSP=TI+(ZERO-A-B)*SK 

Al»=(TIN-TSP)/TSP 

ERR=ABS( A5) 

IF(ERR-XERR)201 ,201 ,203  • 

203  XERR=ERR 

201  CONTINUE 

X0UTA  =  X0UTA  +  2.*3.1 A  1592/180. 
IF ( AO A- 1. ) 299, 304.3  04 

299  IF (XOUTA-3. 1AI592>202,202, 301 

202  A0=(XL2+XL3)*(XL2+XL3) 

AR=XL1*XL1+XL4*XL4-2.*XLI*XL4*C0S<X0UTA) 
IF(A0-AR)301 ,301 ,30  0 


\ 


■ 

'..:•■■     ■-■ 

'  • 

.' 

■■:  ■ 

i       I 

■   ■  •  ■  ' . 

■ 

• 

••'  .'I 


1  304  IF(X0UTA-3. 141592)305,305.401 
j  305  AQ=<XL24XL3)4(XL2+XL3) 

•   '    ..  .■ 

■■■'■-'■     '  ■     . —  ■ 


A«sXLl*XLl*XL4*XL4-2.4XLl*XL4*COSC XOUTA) 

IF(AO-AR)401,401,300 


•  -. 


■    . 
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i 


300."  CONTINUE 

IF( AHA-I . )204,400.400 
204  1F(XEWU-HNO>250.250.252 
250  UNO=XFW«     •  .'■■'•. 

A XL  1= XL  I 

AXL?=XL? 

A  XL  3=  XL  3  _  ^ 

AXL4=XL4 

AZt:«0  =  ZE«0 

AOUTA=/OUTA 

WRITH(3,90  00)AXLl , AXL2 . AXL3 , AXL4 . AZEWO ♦ AOUTA .BNO 

GO  TO  254 

252  WHI  Tl- (3.9005>XL4  ' 
WHITE (3, 0002) XL  I . XL2 « XL 3. XL4 . XERR 
XL4=CXL4+XLINC*C0U 
IF(XL4-i .06*CXL4)302»302.303 

302  GO  TO  100 
301  CONT  JNUE 

253  WR2 TLJ3.9004  )XL4 

254  XL4=CXL4+C0U*XLINC 
IF<  XL4-i.06*CXL4)30  2.302.303 

303  \*RITE{  3.9003) 
Wtf  X TE( 3.900  0 )AXLi . AXL2 .AXL3 .AXL4 .AZERO . ACUTA. BNO 
AUA=1  . 
XL4=AXL4 

cou=o. 

«OUTA=AANG 

GO  TO  !01' 
400  CONTINUE 

404  IF( XFRR-BNQ2)45Q. 450.452 
450  BN02=XEH«   • 

BXL!=XL1 

BXL2=XL2 

BXL3=XL3 

BXL4=XL4 

8ZER0=ZE«O 

BOUTA=ZOUTA 

WWITE(3.9000)BXL1 . BXL2.UXL3 , BXL4,BZER0 .BOUT A ,BN02 

GO    TO    4L>4 

452  WWITfc < 3.9007)Z0UTA 

WHI TEC  3. 9 0 02) XL  1 • XL2 • XL3. XL4 . XERR 
XOUTA=AANG*COU*AlNC 

IF(XOUTA-l  .  01  3*0 ANG)  402.402.  403 
402  GO  TO  101  : 

402  CONTINUE 

453  W«ITE<3,9000)Z0UTA 

454  XOUTA=AANG*COU*AINC 
IF(X0UTA-1 . 013*BANG>4  02.4  02.4  03  '"* 

403  WHITEC3.9009) 
WW[TE(3.9000)BXLi . BXL2.BXL3. BXL4 .BZERO .BOUT A , BN02 
STOP 


.•. 


• 


■. 


■  • 


•.'. 


:    .» 


• 


■ 


I  i 


END 
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SUIWUUT  INE     TOWQ< XL  I , XL2 . XL3 . XL4 . XOUTA, T t N t CSB »CSA ♦ SNAB tCSAB ,£ .F 
1SNA,SNB.XKR.F23X.F23Y»A,U) 
100?     AG=<XL4*XL4+XLI*XL1-2.*XL1*XL4*C0S(X0UTA>) 

17     XL5  =  SGI?T(  AG) 

16    CSNC0=(XL5*XL5-XL3*XL3-XL2*XL2)/«-2.*XL3*XL2) 

AH=l.-CSNCD*CSNCD  '         '  "',     •'.• 

5NCD  =  S0RT(AH)     '  •  '■•"'■■     '  .-'''■''"'   \    \\       '  '  ;      "     ' 

20  SNF=S1N(X0UTAJ*XL1/XL5  •     . 
PP=l. 

21  AA=XL4*XL4  •:•••'    .-,  ;■'     •■:•'_  ■     -  .  .    •  .    '  • 

22  UB=XL5*XL5  ' 

23  CC=XL1*XLI  ' 

24  lF<niUAA-CC>25. 26.27 

25  ^-bNt- 
GU  TO  499  :.;•/•.  ' 

26  F  =  3.  141592/2.  •   i 
GO  TO  28 

27  Z=5NF 


GO  TO  379 


28  5NE=XL2*SNC0/XL5 

■'.'■■.■ 

2Q   DP—  P - 

30  00=XL2*XL2     '  ,1 

31  EE=XL3*XL3  • 

32  IF(DD-EE-B8>33.34.35 

3  3  Z  =  SNE  •  '  •    '•; 

GO  TO  379  •  '  *  • 

34  E=3. 141592/2.  ,  . 

GO  TO  36  .  .   '.: 

•   ..,•■ 

35  2=SNE 

GO  TO  499 

36  XKR=3. 1A1592-X0UTA-E-F  .       -  .    . 
40  TQUT  =  '50.  .   . 

AAAH=1 ./( XL4*C0S( X0UTA)+XL3*C0S( XKH) ) 

AAAI=1 ./(XL4*COS(XOUTA) ) 

AAAK=SIN( XOUTA }/C0S(XgUTA) 

AAAJ=(XL4*SIN(X0UTA)-XL3*S1N(XKR)  >  /  (  XL4*C0S  (  XQUT  A  ) +XL.30C0S  <  XKU  i 

42  F23X=-TOUT*( AAAH-AAAI )/< AAAK-AAAJ) 
AAAD=1  ./(XL<4*SIN<XOUTA)-XL3*SIN«  XKR)  ) 
AAAE=1  ./( XL4*SIN{ XQUTA)  ) 
AAAP=COG(XOUTA)/SIN(XOUTA) 
AAAG=(XLA*C0S(X0UTA)+XL3*C0S(XKR) > • ( XL4* S I N ( XOUT A )-XL3*S I N ( XKR 

43  F?3Y=-T0UT*( AAAO-AAAE)/{ AAAF-AAAG)  .     ' 
SNA=SIN(XOUTA)*XU4/XLS                                 • 
SNB=SNCD*XL3/XH> 

AO=l .-SNA*SNA  • 

AP=1  .-5NB*SNB  ,'  .  ". 

IF<  AA-CC-F1B)44.45.46 

44  CSA-SOMT(AO) 
GO  TO  47' 

45  CSA=0.  " 

GO  TO  47      •  ... 

46  CSA=-SO«T(AO) 

•  !'.•■.  ... 
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47     U0=XL?+XL2 
40     IF(D0  +  UD-E6  )49. 50.51 
49    CSIi  =  -S<iWT<  AP) 
CO     TO    52 

be  csn=c. 

GO  TO  52 
51  CSU=50WT( AP) 
.52^CSAP=CSA*CSH-SNA*SNB  .  :... 
55"SNAH=SNA*CSB+CSA*SNH 

PP=3. 

ir(CSA-0. ) 705. 706. 7  07 

705  Z  =  5NA 
GO  TO  490 

706  A  =  3. 141592/2. 
GO  TO  709 

707  Z  =  SNA  .  .  i  •    ■  •  (  ,' 

708  GO  TO  379 

709  PP=4.  -  '■:'. 
IF(CSB-0. )710.71 1 .712                           * 

710  Z=SNP 

■ 

GO  TO  499 

711  H=J. 141592/2. 
GO  TO  713 


* 

.            ■ 

1     '. 

•'      ' 

■ 

i           ■         •       ■  ■ 

i 

■ 

. 

•      "  •  ■ 

1 

■  • 
■ 

■  ..••■•;. 

'   ; 

-— 

[ 

.  . 

■[    ■;  ,;:.  •  .  !   (,:.:'•     |   ■  ■  ,  '  '.  :;  •  •'■:.  '..  .  :',Y  ■: '.  • 


•  t :.  (  •  ■■  }   i 


712  Z  =  SNR' 

GO  TO  379  j 

379  IF(Z-. 7)380, 380, 401  ' 

380  XN=-1. 

381  XX=0. 

382  XSUM=Z 

3n  3  XSUN"  =  XSUM*(XN+2.)*<XN  +  2.  )  *Z  *Z/  (  {  XN  +  3  •  )  *  (  XN  +  4  .  )  ) 
XX=XX+XSUM 


•   > 


■ 


_  ■- 


398  IF(XSUM-.00001 )400.400,399 


399  XN=XN+2. 


■ 


y.  GO  TO  303 

400  ANGLE=Z+XX 

GO  TO  600  . 


:  . 


401  XN=-I . 
Z=l.-Z*Z 
Z=SOHT(Z) 

xx-o. 

Q        XSU^=Z 

402  Xr,UM=X5UM*(XN  +  2.  )  *  {  XN  +  2  .  )  *Z*Z/  (  < XN+ 3 . ) * ( XN+4 . )  ) 
XX=XX4XSUM 

404  IF<XSUW-. 00001 )405. 405. 402 

405  ANGLt=3. 141592/2. -(Z+XX) 
GO  TO  600  •   , 

499  IF(Z-. 7)500, 500. 520 

500  XN=-S  o 
XX  =  0o 
XSUM=Z 

501  XSUM=XSUM*(XN+2. )*(XN+2. )*Z*Z/( { XN+3 .)*<  XN  +  4. ) ) 
.  XX=XX+XSUM 

'1909  IF<XSUM-.00001)S03,503.50? 


.  i 

. 

1 

• 

■ 

1 

■ 

• 

i 

• 

■  '  . 

■  '  ■  ■ ' 

(XN+4. ) ) 

\ 
i' 

H 

.   ■    . 
• 

■ 

'- 

j 

.     ■  .     . 

•  ... 

. 

■■ 

.  :-«  . 
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502  XN=XN*2. 

CO     TU     'jOI  •  . 

503  ANGLl-=-(Z  +  XX)+3.Ml592 

GU     TO     000  ,  ■       •'  •  '    ■  .' 

520     XN=-l . 


xx  =  o. 


', 


z=i. -z*z 

ZsSQRl (  Z) 
XSUM=Z 
J521     XSUK.=  Xf.UM*(XN-»2.  )*(XN+2«  )*Z*Z/(  (  XN  +  3.  )  *<  XN+4  .  )  ) 
XX=XX+XSUM 


• 


'j?4  1F(  XSUM-. 00001  »525«525»52S 

525  ANGLf =3. 14 1592/2. +  Z  +  XX 

600  IF(PP-1 • )599,60l »603 

60  1  FtANGLE  '■."'..-  : 


fc02  GO  TO  ?«  . 

'603  IFCPP-2. )599.604»606 

599  XOUTA=ANGLE 

604  E=ANGLE 


•  •   .  ■ 


:»  ■-. 


605  GO  TO  36  .  ''.... 

faC6  IF(PP-3. )605»607.609 


. 


607  A^ANGLE 

tOU     GO  TO  709  • 


.' 


• 


- 


. 


609  IFCPP-4. )608. 610.612 

610  U=ANGLE 

■ 

611  GO  TO  713  ... 

612  CONTINUE  . 


• 


713  CONTINUE 

T!N=F23X*XL2*SNAB-F23Y*XL2*CSAB 


.    ;.ji 


.     ■ 


..-■■■ 

RETURN  ' 

ir  pu  r\  *      ■     '  i    ' 


. 


• 


.  .•  ! 


| 

■ 

■ 

.         - 

-    . 


r     -  ■ 


• 


'• 


w 


.  ■ 

V 

■ 


. 

.'•■..■•. 

• 

'       .  • 

. 

'.    • 

■'. 

■     ; 

' 

.   •■ 

• 

- 

.    - 

."'•';  •  •' 


' 


•      Lt\     ' 
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■     •  ■  •  ■  •    i 

■  '  ....  ..... 

..      •  .  ■  ......  •-  :- ••-•        ■.      1 

|  •  ■..,'..  .  ,        . 

-  I  ■  .   •  .     •  • 

y    i  'FUNCTION     ANGLE(Z)  '  •        .     ..' 

J70     IF(Z-. 7)300, 380,401  ■ 

3eo  xn=-i.  .  .  .     .,  ■ 

3ei  xx  =  o.  ;  •  • 

3R2  XSUK=Z  .     . 


383  XSUf/.=  Xf.UM*(XN+2.)*<XN  +  2.)*Z*Z/<CXN  +  3.)*<XN  +  4.)) 

XX=XX  +  XSUM  •        '    '  .• 

398  IF(XSUM-.00001  )400.400,*399 
299    XN=XN+2.  '  | 

GO  TO  383       .   ,  '  .  f  '    •  .   •  C 

400  ANGLE=Z*XX 

...         ' .   . .  ; 

GO  HO  600  •  * 

401  XN=-1  •  .'■  • 


i 
■  . 


7=1  ,~Z*l  • 

Z=f.OftT(Z) 

xx  =  o. 

XSUM=Z 
4  02  XSUM=XSUM*(XN+2. )*(XN*2. )*Z*Z/( <XN*3. )*<XN+4. ) ) 
XX=XX+XSUM 

404  IF(XSUM-. 0.0001  )40S, 405, 402  ,  •• 

405  ANGLE=3. 14 1592/2. -(Z+XX) 
600  ANGLC=ANGLE 

RETUWN 

end  '  •  ;••■  • 


: .  ■ 


* 


i 


FUNCTION  GANGL  (Z) 


499  IF< Z-. 7)500. 500, 520 


500  XN=-1. 


XX  =  0.  ■••'■; 


-• 


- 


XSUM=Z 

501  XSUM=X5UM*(XN+2. )*<XN+2. )*Z*Z/< { XN  +  3 • ) * (  XN  +  4 . ) ) 

XX=XX  +  XSUM  '   .; 

909  iF(XSUM-. 00001 )503. 503, 502 

502  XN=XN  +  2.  j  ■.  . 

go  to  501  • '   ;  .  .'.'i 

|503  ANGLH=-(Z+XX>+3. 141592  ,  ■    . 

GO  TO  600 

520  XN  =  -1  •  '•''. 


•  ■ 


xx  =  o. 

Z=l.-Z*Z 

Z=S0KT(  Z)  ... 

XSUM=Z  .  ' 

521  XSUM=XSUM*(XN+2. > * < XN+2 . ) *Z*Z/ < < XN+3 . > * < XN+4 . ) ) 
XX=XX*XSUM 

524  IFCXSUM-. 00001 >525, 525, 52 i  '  . 

525  ANGLE=3.141592/2.+Z+XX  •   • 

600  ANGLE=ANGLE  .  !. 

■  .'...'-.■■    r  ;. 

WETUMN  v-  ,  ;  . 

END 
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Pages  in  Study:  40 

Major  Field:  Mechanical  Engineering 

Scope  and  Method  of  Study:  The  coordinated  motion  linkage  problem 

can  be  solved  by  an  approximate  graphical  overlay  method.   The 
graphical  method's  results  can  be  then  considerably  refined  by  the— —  >  \ 


application  of  the  variation  of  parameter  technique.  An  example 


: 


problem  which  is  included  in  the  report  considers  a  linkage  which 

must  satisfy  certain  external  torque  requirements.   It  is  analyzed  as 

f 
a  coordinated  motion  problem.  The  linkage's  measure  of  excellence  is 

■ 

the  absolute  value  of  the  difference  of  the  desired  torque  of  link 
two  minus  the  actual  torque  on  link  two  divided  by  the  desired 

torque  on  link  two.  A  computer  program  was  developed  to  perform 

•  i 
the  variation  of  parameter  technique  on  the  graphical  solution.  It 

can  be  easily  modified  to  fit  most  coordinated  motion  problems.       ...... 

•  ,■»■'• 

Findings  and  Conclusions:  The  results  of  this  report  indicate  that  the 
variation  of  parameter  technique  is  a  valid  and  useful  method  for 
reducing  the  error  in  linkage  coordinated  motions.  The  steps  for 
performing  the  technique  are  clearly  outlined.  The  included  computer 
program  is  of  great  use  for  application  of  the  technique.  The  result3 

of  the  example  problem  point  out  that  error  can  be  significantly  reduced  . 

'■-•■••    '•        '■■'.',•/.,,...  i;;;;;)  rxiVl'  OfcTPtf: 

■i  •  "'■'      ■    " 


■(jclicec 


by  use  of  the  technique.  The  computer  was  utilized  to  reduce  time 
and  increase  accuracy  of  the  calculations.  The  included  program  can 
be  modified  to  make  an  even  better  approximation  for  the  calculation 
of  the  maximum  error. 
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